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Abstract 

Most  problems  involving  dielectrics  cannot  be  solved  exactly,  e.g«, 
diffraction  by  a  dielectric  wedge.  It  is  the  purpose  of  this  work  to 
describe  a  technique  by  which  the  transition  conditions  at  an  interface 
between  tvzo  dielectrics  can  be  conveniently  and  accurately  replaced  by  a 
single  mixed  boundary  condition.  The  method  is  useful  when  one  is  interested 
in  the  field  in  only  one  of  the  dielectrics*  The  boundary  condition  takes 
the  form 

|H  +  Au  +  B  ^  =  0 
^^        as^ 

where  v  and  s  are  the  surface  normal  and  tangent  respectively.  When  a 
single  interface  is  present,  so  that  there  is  no  diffraction,  it  is  proved 
that  the  maximiun  percentage  error,  over  the  entire  far  field  of  an  arbitrarily 
placets  line  source  is  small.  For  example,  we  can  tabulate  the  maximum  per- 
centage error  for  transverse  electric  excitation  as 


Index  of 
refraction  =  n 
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The  percentage  error  becomes  very  large  for  n -%-  1  because  for  n  =  1 
there  is  no  reflection  at  the  interface;  however  the  absolute  error  remains 
small  and  bounded. 

Additional  confidence  in  the  technique  is  acquired  by  using  the  boundary 
condition  to  replace  a  dielectric  interface  in  a  non-trivial  problem  which 
does  involve  diffraction.  A  comparison  of  the  results  of  the  approximate 
version  and  the  exactly  formulated  problem  shows  that  there  is  excellent 
agreement  between  the  two. 
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Introduction 

Most  electromagnetic  problems  involving  more  than  one  dielectric 
Medium  are  not  capable  of  exact  solution.  An  example  of  such  a  problem 
is  sketched  in  Figure  1,  One  is  interested  in  obtaining  a  solution  of 
the  wave  equations 

(V^  +  k^)  u(x,y)  =0,  i  =  1,  2,  3 

for  the  indicated  geometry  when  a  plane  wave  is  incident  in  medium  I. 
It  is  assumed  that  the  area  of  prime  interest  is  region  I,     The  boundary 
conditions  to  be  satisfied  are  continuity  conditions  determined  from 
Maxwell's  equations.     As  the  problem  stands,   it  is  not  solvable  by  available 
mathematical  techniques.     However,  one  notes  that  the  two-part 
boundary-value  problem  of  Figure  2  is  amenable  to  the  Wiener-Hopf 
procedure  if  BC, „  and  BC_     are  mixed  boundary  conditions  with  constant 
coefficients. 

The  following  procedure  then  suggests  itself  for  obtaining  an 
approximate  field  in  region  I: 

(1)  Replace  the  dielectric  interface  between  regions  I  and 
II  by  a  linear  boundary  condition  BC--  =  C,   and  do  the  same  for  the 
pair  of  regions  I  and  III. 

(2)  Solve  the  problem  of  Figure  2  exactly  to  obtain  an 
approximate  solution  to  the  problem  of  Figure  1  in  region  I, 

The  success  of  this  approach  depends  upon  constructing  a  boundary 
condition  to  replace  a  dielectric  interface  and  validating  the  procedure. 


il 


In  the  present  report  (Part  I  of  two  parts)  we  discuss   the  formulation 
and  justification  of  the  new  boundary  condition,   and  in  a  subsequent 
report   (Part  II)  we  solve  the  problan  of  Figure  1  as  an  illustration 
of  its  use. 
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Figure  1:     Idealized  Three-Media  Problem 
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Figure   2:   Reformulated  Approximate   Version  of  the   Problem 
Shown  in  Figure   1   in  Region  I. 


In  section  1  we  determine  the  coefficients  in  the  new  boundary 
condition  by  comparing  the  reflection  coefficient  it  defines  to  the  actual 
Fresnel  reflection  coefficient  appropriate  to  a  dielectric  surface.    W© 
present  numerical  data  that  shows  that  there  is  little  error  in  approximating 
the  Fresnel  coefficient  if  our  methods  are  eii?)loyed.     Certain  limitations 
and  remarks  concerning  this  procedure  are  presented  in  Section  2. 

Sections  3  and  h  are  devoted  to  an  investigation  of  the  utilization 
of  the  new  boundary  condition  at  the  interface  of  a  dielectric  half-space. 
In  these  sections  it  is  shown  that  the  introduction  of  the  boundary 
condition  leads  to  a  problem  with  a  unique  solution  vrtiich  satisfies 


±r 


reciprocity.     Once  reciprocity,  uniqueness  and  existence  have  been 
demonstrated  we  are  in  a  position  to  assert  that  this  boundary  condition 
leads  to  (l)  well-defined  problems,  and  (2)   guarantees  that  the  eiror  in 
the  far-field  of  the  approximate  vei^ion  of  the  standard  problem  of  the 
radiating  line  source  above  a  dielectric  interface  is  small. 

Since  the  problem  treated  in  sections  3  and  k  involves  no  diffraction, 
confidence  in  the  utilization  of  the  new  boundary  condition  is  amplified  in 
Section  5,   in  which  we  compare  an  exact  solution  of  a  diffraction  problem 
involving  two  dielectrics  and  its  approximate  counterpart.     We  find 
strikingly  good  agreement  between  the  two  results* 
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1,     Construction  of  the  Boxindaiy  Condition 

In  order  to  replace  a  dielectric  interface  by  a  linear  boundary 
condition  it  is  necessary  to  determine  what  characterizes  the  interface 
so  far  as  wave  propagation  is  concerned*     If  one  solves  the  two-medium 
problem  of  a  line  source  above  a  dielectric  interface  (Fig.  3)  then  one 
obtains  integrals  of  the  form  (cf.  Bremmer,    [8]). 


(1) 


j    R(e)  e^ 


ilq^cos  6 


de 


in  addition  to  the  field  that  would  exist  in  the  absence  of  the  boundary. 
In  the  integral  of  (l),  R(&)  is  the  appropriate  Fresnel  reflection 
coefficient  and  the  contour  C  is  the  familiar  contour  defining  the 
Hankel  function. 
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Figure  3»  Contour  defining  the  Hankel  function 
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The  diffracted  field  represented  by  (l)  then  has  the  interpretation  of 
a  summation  of  plane  waves  traveling  in  all  directions,  real  and 
imaginary,  which  have  as  a  weight  factor  the  Fresnel  coefficient  defined 
in  the  complex  9-plane. 

One  can  also  observe  that  a  linear  boundary  condition  at  the 
interface  will  define  a  reflection  coefficient  if  we  seek  plane  wave 
solutions.  When  medium  (2)  is  almost  conducting,  but  only  then,  it  is, 
in  fact,  known  that  a  boundary  condition  of  the  form'- 

u  +  \u  ■  0 

y 

represents  the  phenomena  fairly  well.     Therefore,  for  a  more  general 

case  we   seek  to  include  more  parameters  so  as  to  improve  the  approximation. 

For  example,   notice  that  the  boundary  condition 

(='       (i^  *  *  *  ^  ^'"-°'       ^•° 

will  generate  for  the  wave  equation 


(3)  (V^  +  k^)  u(x,y)  =  0  ,  y  >  0 

the  reflection  coefficient 

o 

cos  »  -  (A-  B  sin    &) 

(U)  R(e) 5— - 

cos  9  +  (a-  B  sin    9) 

where  &  is  the  angle  of  incidence  of  a  plane  wave  (cf.  Fig.  i^  ) 


-  3  - 


"inc" 


ik(x  sin  9  -  y  cos  B) 


incident  upon  the  boundary  at  y  ■  0  . 
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Figure  U»  Plane  wave  incident  and  reflected  at  y  -  0. 


To  derive  the  expression  (U)  one  assumes  a  geometricsally  reflected  field 


ti     -  -  Re 

ref 


ik(x  sin  &  +  y  cos  &) 


and  the    amplitude  coefficient  R(e)  is  determined  so  as  to  make  the  total 
field 


^inc  *  ^  Vf 


satisfy  the  boundary  condition  (2). 
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An  approximation  procedure  is  then  evident:  Find  a  boundary 
condition  whose  reflection  coefficient  is  a  suitable  approximation  to  the 
appropriate  Fresnel  reflection  coefficient.  One  is  reminded  that  if  u 
represents  E_  the  Fresnel  coefficient  is  ( [$] ,   p.  U92  ff . ) 


(5)R^- 


cos  e  -  .jly  l-(ri]  sin^  § 


^|V^ 


^i"(^^«)  k_/   ,^2  ,  2 


cos  ©  ♦  rr=' W  1-  (  p^  1  sin  6 


vrfiere  ©     is  the  angle  of  refraction  and  k^  is  the  propagation  constant 
in  the  lower  media  y  <  0.     If  u  ■=  H     then  the  reflection  coefficient  is 


■^:V^^  -' 


In  (5)  and  (6)  we  assume  the  magnetic  peimeabilities  li^i^  and  n^  to  be 
identical  for  convenience  only.  If  we  conpare  (5)  and  (6)  with  (U)  then 
we  see  that  we  can  approximate  (5)  and  (6)  at  least  for  all  real  angles 
of  incidence  by  choosing  A  and  B  properly.  To  make  this  selection  we  note 
that  if 


(7) 
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then  the  approxijnation  (U)  will  agree  exactly  with  either  (?)  or  (6) 
for  nonnal  incidence  (i.e.,  9  =  O).     For  glancing  incidence  (©  =  n/2), 
the  expression  (U)  will  agree  with  (5)  or  (6)  regardless  of  the  choice 
of  B  since  both  are  equal  to  -1  for  6  »  n/2. 

A.  Angular  Matching 

We  can  choose  B  in  many  ways  to  improve  the  approximation.     For 

example  B  may  be  chosen  so  as  to  make  the  approximation  exact  for  some 

angle  9,  where  0  <  6.  <  n/2.     There  is  nothing  to  suggest  that  any  angle 

S-   is  to  be  preferredj  for  the  balance  of  Section  1  we  shall  discuss 

possible  choices  and  their  motivations.     For  any  choice  of  9,   the 

Bmn.  n»#  ^^6  determined  by  matching 
TEjiM 


B„ 


1  .  "'^;™  sin^      with  -/(l-(l/n)2sin^. 


where  n  =  k^/k.  is  the  index  of  refraction  or 


(8) 


^TE,TM^®1^ 


^TE, 


TM 


1  -/(l-(l/n)^3in\) 


sin^e. 


B.  Binomial  Matching 

If  n  >  1  we  can  expand  the  radical  in  (8)  by  the  binomial  theorem 
to  obtain  from  (7) 


(9) 


TE 


^TM 
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Note  from  (9)  that  the  leading  terms  of  B^t-  and  B_,  are  independent  of 

ill      wi 

e  ,  the  fitting  angle.  For  analytical  purposes  it  is  very  convenient  to 
choose  the  B's  as  just  these  leading  terms 


(10) 


TE  "   2n 


'»■  ^ 


Ihs  selection  (10)  will  be  referred  to  as  binondal  matching.     Physically 
it  corresponds  to  improving  the  match  of  the   reflection  coefficients  in 
the  neighborhood  of  6  =>  0,  or  normal  incidence.     This  is  illustrated 
in  Graph  1  where  we  plot  the  percentage  error  between  the  reflection 
coefficients  (5)  and  (6)  and  their  approximations  (U)  with  binomial -match 
coefficients  for  n  »  1.6.    Note  that  in  the  neighborhood  of  Brewster's 
angle  9g  ,  R(Sg)  "  0  and,  hence,  the  percentage  error  (but  not  the 
absolute)  for  TM  approximation  is  necessarily  unbounded.     In  Graph  2 
we  plot  the  absolute  errors  and  we  see  that  the  absolute  errors  in  the 
vicinity  of  B-  remain  small.     Indeed  the  absolute  errors  for  both  the  TE 
and  IM  approximations  have  about  the  same  values, 

C,  Brewster's  Angle  Matching 

We  can  eliminate  the  unbounded  percentage  errors  for  the  TM  case 
In  the  vicinity  of  Brewster's  angle  by  choosing  B  so  as  to  yield  a  perfect 
match  at  Brewster's  angle  9«  ■  tan"  n  in  which  case  (8)  becomes 
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Curves  3  and  h  illustrate  the  percentage  and  absolute  errors  for  this 
choice  of  ©^  for  n  »  1,6. 


D.  Chebyshev  Matchings 

Still  another  cl»ice  for  B  is  possible  at  least  for  the  TE 
approximation J  this  selection  is  a  value  of  B  so  chosen  as  to  minimise 
the  maximum  percentage  error  for  the  TE  reflection  coefficient.     We  note 
that  the  magnitude  of  the  error  changes  sign  for  some  value  9^| 


if  we  choose  a  value  of  9»  such  that  the  raaxiimim  positive  percentage  error 
equals  the  maxLmura  negative  percentage  error  then  we  shall  have  the  least 
raaxiinum  error  for    the  intenral  0  <  &  <  ii/2  by  the  ChebysheV  criterion. 
The  mathsBiatical  analysis  to  find  the  optimum  Chebyshev  matching  would  be 
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interesting.     Here  we  have  done  it  by  a  trial  and  error  processj   the 
values  of  d_  required  to  afford  this  match  for  the  TE  case  are  plotted  in 
Graph  5,     The  maxiBiuin  associated  errors  for  the  reflection  coefficient  are 
plotted  in  Graph  7  as  a  function  of  n.     It  is  not  possible   to  repeat 
this  procedure  for  the  1M  case  owing  to  the  unbounded  nature  of  the 
percentage  eiror  in  the  vicinity  of  Brewster's  angle,  however  a  Chebyshev 
matching  is  possible  for  both  the  TE  and  TM  cases  by  minimizing  the 
absolute  error,  but  we  have  not  carried  this  out. 

It  seams  reasonable  to  use  the  boundary  condition  (U)  with 
Brewster's  angle  coefficients  for  problems  of  transverse  magnetic  excitation. 
For  a  problem  of  transverse  electric  excitation  we  can  use  the  values  of 
©_  in  curre  (5)  in  formula  (8)  to  yield  a  Chebyshev  fit.     The  maximum 
percentage  errors  that  arise  for  real  angles  of  incidence  for  these  two 
approximations  are  plotted  for  reference  in  cuirres  6  and  ?  as  a  function 
of  n,  the  index  of  refraction. 

Even  better  Chebyshev  fittings  are  possible  by  relaxing  the  demand 
of  a  perfect  fit  at  &  =  Oj  this  permits  choosing  both  A  and  B  to  improve 
the  match. 
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This  procedure  would  yield  a  perfect  fit  for  four  values  of  6  over 
the  interval  0  <  9  <  n.     This  particular  matching  was  not  carried  out, 
owing  to  the  labor  involved  in  doing  it  empirically,   and  the  excellent  fit 
of  the  cruder  approximations. 

E.  Grazing  Matching 

For  certain  problems  an  alternate  procedure  is  available  if  a 
physical  inquiry  is  of  such  a  nature  that  most  of  the  field  of  interest 
is  associated  with  a  particular  angle  of  incidence.     An  example  is  a 
traismitter  and  distant  receiver — both  located  near  the  earth's   surface. 


The  geometry  of  this  problem  is  such  that  one  would  expect  most  of  the 
energy  arriving  at   the  receiver  to  be  associated  with  the  arrival  of  rays 
in  the  neighborhood  of  grazing  incidence.     For  this  problem  one  wo\ild 
choose  a  value  of  B  to  match  the  derivative  of  the  reflection  coefficient 
at  grazing  incidence.     From  (8)  we  see  that  the  requisite  B(n/2)  would  be 


^^^  ^T^,™  ^"/2)  -  A^^^ 


-v^ 
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In  subsections  A-E  we  have  indicated  several  possible  choices  of 
the  coefficients  in  the  B.C.  For  the  balance  of  the  report  we  shall 
assvune  that  the  coefficients  have  been  chosen  for  a  particular  problem,  and 
confine  our  discussion  to  the  question  of  what  is  true  for  any  such  choice 
of  A  and  B^ unless  otherwise  indicated. 

F.  Genera!! izations 

Beside  the  above-mentioned  methods,  other  refinements  are  possible. 
Thus  any  desired  degree  of  accuracy  can  be  obtained  by  adding  additional 
terms  to  the  boundary  condition:     By  adding  such  terms  as 

9  u  8  u 

8x  3x 

to  the  boundary  condition  (U)  we  would  have  generated 

sin%  ,         sin  9,  ... 

in  the  approximate  reflection  coefficient.     The  coefficients  of  these  terms 
would  then  be  available,  either  to  match  as  many  terms  in  a  series 
expansion  of 

Vl  -  i-f  sin^  9 

"  n 

as  needed  for  superlative  accuracy,  or  for  use  in  a  Chebeyshev  approximation. 

2,  Limitations  and  Remarks 

Our  approximation  procedure  is  limited  by  the  restriction 
k_  >  k,,  or  n  >  1,     For  if  k^  <  k,   then,   there  arises  a  real   angle  of 
incidence  corresponding  to  total  reflection  for  irtiich  the  second  radical 


-  n  - 


in  (1.5)  and  (1.6)  becomes  imaginary.     Since  (l.U)  involves  constant 
coefficients  the  boundary  condition  is  unable  to  reproduce  the  transition 
from  a  real  reflection  coefficient  to  a  complex  one  for  real  angles  of 
incidence.     This  Is  clear  by  inspection  of  the  reflection  coefficient 
(l.U).     So  long  as  both  A  and  B  are  .isaiy   the  reflection  coefficient 
remains  real  for  any  real  angle  ©•     If  A  and  B  are  imaginary  or  complex, 
then  R(8)  will  always  be  complex  for  any  range  of  0. 

One  can  also  observe  fran  inspection  of  the  curves  6  and  7  that  the 
approximation  deteriorates  for  n  ^  1.     Physically,  this  corresponds  to  an 
interface  with  a  reflection  coefficient  in  the  vicinity  of  zero  and  hence 

difficult  to  approximate  on  a  percentage  basis.     However,  this  fact  allows 

[91 
one  to  use  other  approximate  methods'-  -'.  However,   the  absolute  errors 

remain  small  and  bounded.     As  a  consequence  the  boundary  condition  (l.U) 

specialized  for  n  ■  1  may  be  useful  as  a  nodel  of  a  perfect  absorber  since 

the  reflection  coefficient  it  defines  is  so  small  for  most  of  the  range. 

A  curve  of  this  case  is  drawn  in  Graph  8.     It  is  unavoidable  that  the 

reflection  coefficient  goes  to  -1  for  ©  »  n/2.     This  behavior  is  forced  by 

the  form  of  (l.U)  which  is  true  even  in  the  exact  case  if  n  /  1. 

A  special  feature  of  the  new  boundary  condition  is  that  it  involves  a 
second  derivative  in  x.     In  Section  h  we  show  that  introduction  of  svich  a 
term  does  not  violate  reciprocity,  and  the  method  of  demonstration  will 
show  that  we  can  include  additional  derivatives  with  respect  to  x  of  even 
order  without  sacrificing  reciprocity. 
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3.   Existence  and  Uniqueness 
A.  Existence 


In  this  section  we  shall  obtain  an  explicit  solution  to  the 
following  problem 

(1)       (V^  +  k^)  u(x,y)  -   Ui6(x)  6(y-h),        y  >  0 


p  =  y(x  +  y  )   J       ©  =  arctan  (y/x) 
Once  a  solution  to  this  problem  has  been  obtained  and  discussed  we  shall 

show  that  we  have  the  only  such  solution.     The  purpose  of  this  demonstration 

and  that  of  Section  h  in  which  reciprocity  is  proved  is  to'  show  that  the 

introduction  of  the  new  boundary  condition  (2)  leads  to  a  well-defined 

approximate  formulation  of  the  standard  problem  of  the  radiating  line 

source  above  a  dielectric  interface. 

A  solution  to  the  above  problem   can  be  obtained  very  simply  by 
expressing  the  excitation  field  u  (x,y)  by  separation  of  variables,   in 

G 

the  form 

n\  -^    f     ivx  +  ±)/k^  |y  -  h| 

(h)    u^(x,y)  =  ^y-h^pj  •=  -  i        5 — = dV 

e     '  o  h  n  )  i-x w 


irtiere  p.    =  yx    +  (y-h)     and  expressing  the  secondary  image  field  u  (x,y) 
in  a  similar  form,   namely, 


page  22. 
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.iv^?T:3 


Here  G(v)  is  to  be  determined  by  the  condition  that  the  total  field  (6) 

(6)   u^(x,y)  -  u^(x,y)  ♦  Ug(x,y) 

satisfies  the  boundary  condition  at  y  ■  0.  The  contour  C  is  taken  in  the 
complex  v-plane  along  the  real  axis  as  shown  in  Fig,  8  which  also  indicates 
the  choice  of  branch  cuts  for  Vk  -  v  .  In  this  section  we  assume  k  to 
be  real  unless  otherwise  mentioned. 

If  we  apply  the  boundary  condition  (2)  to  the  combination  (6)  we 
find  that  if  we  choose 


(7)      G(v) 


^ 

■7 

-  (A 

2 

^F 

v2 

♦  (A 

then  the  boundary  condition  (2)  will  be  satisfied.     It  will  follow  frcro 

the  asynptotic  development  of  the  integrals  (5)  that  the  radiation  condition 

(3)  will  be  satisfied.     The  complete  solution  can  be  displayed  as 


(8)  u^(x,y)-H^^^(kPj^)  -i       /    ■■■     ,         , ^  •  2 —^^ dv  . 


k 

■7 

(A 

Bv^ 
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Figure  ^:  The  contour  C  in  the  v-plane  defining  (8).     Tn  this 
illustration  the  imaginary  part  of  k  is  neglected. 
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In  order  to  discuss  the  integral  (5)  that  appears  in  (8)  it  is  convenient 
to  introduce  the  polar  coordinates 


(9a) 


P_jj  -   7  x^  +  (y  +  h)^ 


a      -    tan' 


-1    2_^_h 


or 


(%) 


X  -  p  .    cos  © 


y  +  h  ■  p  ,    sin  6 


It  is  also  convenient  to  make  a  transformation  from  the  v-plane  to  the 
cos  (^  plane;  where 


(10) 


V  ■  k  cos  fi  , 


When  this  is  done,  then,  in  the  region  x  >  0,  y  >  0,  the  integral  (5) 
becomes 


(11)    Ug(p_j^,  e)  -  i   I   R(0)  e 

C'(0) 


ikp^j^cos  (0-©) 


vhere  R(0)  is  the  approximate  reflection  coefficient  (1*U)  extended  to  the 
complex  0-  plane 


(12)     R(^)  -  3in  ^  -  (A  -  B  cos^)   ^ 
sin  0  ♦  (A  -  B  cosl^) 
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Owing  to  the  restricted  choices  of  A  and  B  as  discussed  in  this  report 
(Section  1),  there  are  no  roots  of  the  denominator^ of  R()d)  for  real  0  in 
the  range  0  <  0  <  n/2.     Consequently  we  can  apply  the  method  of  saddle 
points  (see  Part  II  for  a  complete  discussion)  directly  to  (ll)  to  yield 
the  asymptotic  development 


(13)         n^ip,^,  B)  -Vl^    ^^^^  « 


'/ 


"=<^-^-l'.o 


"'^-h'  ' 


From  the  development  of  the  Hankel  function 


(li.)   B^'Wt)-|^.'^''^"^^0 


(kp,) 


it  follows  that  the  entire  field  at  large  distances  is 


(15)    u(p,e)  'J^      1  +  R(©) 


"^^P-h-Ph^ 


i^ip^^p 


plus  terras  of  order  p.      .     For  p  .    »  1  we  have  that 


(16)       p^^  -  pj^    -    h  sin  e 


so  that  (15)  becomes 

(17)      u(p,&)  'J-^    [l  +  R(») 


ikh  sin  © 


e 


lk(pj^  -  n/k) 


which  we  see  is  just  (U.5)  within  a  constant  factor  if  we  note  that  the 
angle  used  in  this  section  is  the  complement  of  the  angle  used  in  Section  U. 


see  Appendix  C 
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Hence  the  field  so  obtained  possesses  the  predicted  properties  of 
equation  (3)  of  this  section.     Equations  (1),    (2)  have  been  fulfilled 
automatically  in  the  course  of  the  derivation. 
B.  Unigneneas 

We  can  utilize  the  far-field  behavior  (1?)  of  u(x,y)  to  prove 

the  uniqueness  theorem:  Let  u(x,y)  be  any  function  which  satisfies 


(18)        (7^  +  k^)  u(x,y)  -    6(x-Xq)  eCy-y^)  , 


y  >  0 


a(x,o)  -  0  ,  y  ■  0 


where  A  and  B  are  arbitrary  real  constants. 

In  addition  we  shall  stipulate  that  u(x,y)  can  be  developed  as 


g(©) 


jikp 


+  0 


(kp) 


3 


(20)  u(p,9)-< 


0(1) 


deM 


de!f 


for  kp  »  1  where  M  is  some  set  whose  coirplement  M  la  a  set  of  measure 
zero  on  the  interval  0  <  9  <  n  . 

Proof 

Let  u(x,y)  =  u,(x,y)  -  Vi-iy.fj)   where  u.  and  vl^   *^  ^^  **°  functions 

that  satisfy  (18),  (19)  and  (20).  *rhen  u(x,y)  is  a  solution  of  the 

homogeneous  wave  equation 


(21)    (V^  ♦  k^)  u(x,y)  -  0 
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and  in  addition  satisfies   (19)  and  (20).     Construct  the  function  T(x,y) 
in  terms  of  u(x,y)  as 


(22) 


r(x,y)  gl^  +  Au  +  B^,  y>0 

^    "^  ax 


then  v(x,y)  satisfies  Dirichlet  conditions  at  y  »  0  as  well  as  the 
honwgeneous  wave  equation  for  y  >  Oj  we  shall  first  show  that  v  must 
vanish  identically  for  y  >  0,     This  can  be  seen  very  simply  by  an 
application  of  Green's  theorem  about  a  semi-circle    f""  o^  radius  p  to 
T(x,y)  and  its  conjugate  v  (x,y). 


Namely  we  have 


(v*  g  -  V  |l!)  ds  -   f      {y*V^y  -rvV)  dir    , 

A(r) 


and  by  virtue  of  the  boundary  condition  (19)  (which  implies  that  v  (x,o) 
v(x,o)  5  0)  it  follows  that 
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(23)  (v*  1^  -y  1^)  pde  -  J        (vVt  -  ry^v*)  dt 

-b  ACT) 


2 

If  k    has  a  positive  imaginary  part  & ,  the  line  integral  will  decay 

exponentially  with  increasing  p  leaving  us  with 
lim  2ior'     J       v^v  dt    =  0 

p->00  /(pj 


2 

which  implies  v  b  0,     If,  on  the  other  hand,  k     is  real,  the  integral  over 

A(P)  vanishes  identically  which  leaves  us  with 


Once  again  if  we  let  p  ->oo  we  can  use  the  condition  (20)  to  change  (2li) 
to  the  form 

(25)       J     p(©)  p*(&)  de  =  0 

where  p(&)  is  the  pattern  function  for  v(x,y)i  this  in  turn  implies 
p(&)  =  0  a.e.     The  function  v(x,y)  therefore  has  a  vanishing  far-field 
amplitude  in  its  region  of  definition  y>  0.     We  wish  to  apply  a  theorem 
of  Rellich  in  order  to  show  v  s  0.     The  relevant  theorem  applies  to  a 
function  regular  outside  some  circle.     We  therefore  proceed  to  extend  v 
as  follows:     Since  v(x,y)  vanishes  on  the  x-axis  it  can  be  analytically 
continued  into  the  lower  half-plane  by  the  reflection  principle.     Hence 
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V  is  now  regular  in  the  complete  x,  y-plane.     We  are  now  in  a  position 
to  apply  Rellich's  Theorem.     This  theorem'-  -'   states  that  zero  is  the 
only  wave  function  which  is  regular  in  the  entire  x,y-plane,   and  which 
vanishes  more  rapidly  than  0(— )   for  p  -^oo  and  ©cl  where  I  is  any  interval, 
Equation  (25)  therefore  implies  thn-.  ^    3  0. 

These  results  imply  that  any    solution  u(x,y)  of  the  homogeneous 
wave  equation  with  the  boundary  condition  (19)  satisfies  the  differential 
equation. 

(26)         1|H+Au  +  B^  =  0 

^  ^y  ax^ 

for  all  y  >  0,  It  is  a  simple  matter  to  show  that  this  is  inconsistent 
with  the  assumed  aymptotic  development  (20)  thus  yielding  the  desired 
contradiction.  To  obtain  this  result,  substitute  (20)  into  (26)  and  we 
obtain  almost  everywhere 


ikp 
(27)   g(e)  5 


k  sin  9  +  A  -  B  cos  9 


in  which  x  ■  p  cos  &  and  y  =  p  sin  9.  Equation  (2?)  can  be  satisfied  only 
if  g(9)  vanishes,  or  else,  if 


(28)   k  sin  &  +  A  -  B  oos^  »  0, 


a.e« 


Since  (28)  is  a  quadratic  in  sin  8  it  follows  that  the  equation  (28)  can 
be  satisfied  for  at  most  two  values  of  sin  9.  Hence 


(29)    g(9)  =  0    a.e.,    0  <  9  <  n. 
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The  result  (29)  is  not  sufficient  to  imply  that  u(x,y)   =  0.     We  shall 
complete   the  theorem  by  an  appeal  to  Stokes's  theorem.     Since  u(x,y) 
is  a  wave  function  the  integral 


(30) 


(uV^u*  -  ti*7^u)  dli     -     0 


A(r) 


vanishes  for  any  region  bounded  by  a  curve   f""    in  "the  upper  half  plane. 
In  particular  if   H    is  a  semicircle  of  radivis  p  we  can  compute  the 
integral   (30)  in  another  way.     We  observe  that  if  the  Laplacian  operates 
on  any  wave  function  satisfying  (26)  we  obtain 


„2  ,2  1    Su  ^  B  d  u 

V  u  -  -k  u  -  ^  3^  *  J  -^ 


so  that  (30)  becomes 


A(r)L 


Atn 


-  ^  (U  Uy  *  U  Uy)   *  ^  (uu^  -  '^xx  ^   ^ 


df 


3    3/*  »x 


1    8    /     *x 


dlK 


wheiiOe,  by  Stokes's  theorem 


r 


J 


J  (uu^  -  u^u  )  dy    +    -jj-  uu    dx 


or 


0  -  i^    .    I2    *    I 


-  22  - 


where 


(uu    -  u  u  )  p  COS  8  de 

X  X 


/ti 


T  1  » 

Ig  °  -  -rr  I     11^    P  sin  9  d& 


-P 


As  p  -♦oo  the  first  two  integrals  L.  and  !«  vanish  since  by  (29)  both 


u  and  u     behave  like 
7 


3 
p       ,  a.e.         as   |p|  ->oo        0  <  6  <  n   . 

We  are  left  with 

r 

J        uu     dx     -      CD 
^00 

30  that  u  must  vanish  on  the  boundary,  y  »  0.     Hence  u(x,y)  s  0,  y  >  0 
by  the  same  argument  previously  applied  to  the  function  v. 

Ii,    Validity  of  the  Approximation  (with  respect  to  the  field  of  a  line 
source ), 

In  Sections  1  and  2  we  have  given  no  justification  for  the  boundary  con- 
dition (1.1;)  other  than  it  guarantees  that  the  error  in  the  reflection  coeffi- 
cient is  small  for  plane  waves  at  real  angles  of  incidence.     The  purpose  of 
this  section  is  to  prove  that  if  we  have  a  line  source  above  a  dielectric 


-  22a  - 


R(r^,ry) 


Figure  6:     The  Geometry  of  a 
source  at  S,  a  receiver  at  R 
and  its  imaf^e  at  R*. 
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interface  and  we  replace  the  interface  by  the  boiuidary  condition  (l.U)  then 
the  maxioimn  percentage  error  that  can  arlee  in  the  far-field  by  this 
replacement  is  8iii?)ly  related  to  the  maximum  error  in  approximating  the 
reflection  coefficient.     The  proof  is  based  upon  reciprocity:  If  one 
has  a  source  S  at  (0,h)  above  a  dielectric  interface  at  y  •■  0  and  a 
receiver  atR"(r,r)in  the  far- field  then  the  field  at  R  due  to  the 
source  S  is  of  the  form 


(1)        Ug(S)  -  g(e)  2 +    o(p--^/2j 


and  since  reciprocity  applies  we  can  interchange  the  source  and  i^ceiver 

ikp  _ /rt 

(2)  u-(R)  -  g(e)  2 ♦  o(p-3^2) 

with  the  source  in  the  far-field.     For  a  distant  source  we  know  that  the 
received  field  behaves  like  a  plane  wave  incident  at  an  angle  6  plus 
R(&)  times  a  geometrically  reflected  wave  where  R(0)  is  the  Fresnol 
reflection  coefficient,  i.e. 

ikp  ikp  ,/. 

(3)  ii.(R)  -  2 —    ♦    R(e)  S —   +    o(p-5/2) 

s         y-p"  VT^ 

t  t 

in  which  p    is  the  image  distance  from  3  to  R  ,     Hence  by  reciprocity 


(Ug(R)  »  Ujj(3)     or 


(li)     Up(3)  -  2 +     R(e)  2 +  0  (p-3/2  ) 

vp  vy 


-2U  - 


t 

For  p  »  1  we  can  neglect  the  difference  p  -  p  in  the  radical  but  must 


preserve  it  in  the  exponential 


p  \:^  p  ♦  h  cos  &  ,  p  »  1 


so  that  (li)  can  be  rewritten 

.ikp 


1  +    R(©)  e^^  ^^  * 


(?)      u^(S)  . 

and  if  we  compare  (5)  with  (l)  we  see  that 


(6)    g(e)  -  1  >  R(e)  e^  °°»  « 


Assuine  for  the  moment  that  the  boundary  condition 
(l.U)  satisfies  reciprocity.  Then  the  same  argument  would  yield 

(7)    g'(&)-l*R'(e)e^°°«^ 

as  the  angular  variation  of  the  far-field  if  we  replace  the  dielectric 
interface  at  y  =  0  by  the  boundary  condition  (l.U).     In  (?)  R  (©)  is  the 
approximate  reflection  coefficient  defined  by  (l.U).     Consequently  if 
the  boundary  condition  (1.9)  satisfies  reciprocity  we  can  compare  (6)  and 
(7)  and  make  the  assertion:     Replacement  of  a  dielectric  interface  by 
its  approximation  (l.U)  creates  an  error,  in  the  far-field  of  a  line 
source, which  is  proportional  to  the  error  in  matching  the  reflection 
coefficient  for  plane  waves  at  real  angles  of  incidence. 

We  shall  prove  the  required  reciprocity  theorem  which  validates  the 
above  argument. 
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Proof  of  Reciprocity 

Define  two  solutions  of  the  ware  eqaation  tu  and  Ug 

(8)  (v2  ♦  k2)t.j^  -  6(R) 


(9)  (y^  ♦  k\  -  6(S) 


which  both  satisfy  the  general  boundaiy  condition  at  y  ■  0 

(10)      (-^   *   «   *    P^^    '■■° 


and  the  radiation  condition  at  infinity 


corresponding  to  a  suppressed  time  factor  e~   ,  R  and  S  are  diorthand 

for  the  coordinates  of  the  receiver  (r  ,r  )  and  source  (a  ,s  ) 

X*  y  X*  y 

6(R)  -  6(x  -  r^)  6(y  -  r  ) 
X  y 

6(S')  -  6(x  -  8^)   5(y  -  3   ) 

*  y 


Apply  Green's  theorem  to  u^  and  Up 


W^S-V'V^^-  jK^  -^5?)<^^ 


uhez«  y  is  the  area  bounded  by  the  large  semi-circle  and  S  is  its  boundary* 
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»    X 


As  the  radius  p  of  the  semi-circle  goes  to  infinity  we  obtain  (Appendix  C) 


,00 


2 


Uj^(S)  -  Ug(R)  =  -p    I       (u^ ^  -  u. 


=\ 


)  dx 


-00 


ax 


ax 


since  the  integration  over  the  arc  vanishes  by  virtue  of  the  radiation 
condition.     The  remaining  integrand  is  a  perfect  differential 

alE  W  -^  -  ^  -ai^ 


so  this  integral  vanishes  too  by  the  radiation  condition  thus  proving 
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reciprocity  i.e« 

tljj(S)       -     tlg(R)  . 

As  shovn  in  the  previous  section,  one  is  then  justified  in 
stating  that  approximation  of  the  reflection  coefficient  for  real 
angles  of  incidence  corresponds  to  approximating  the  far-field  of  a  line 
source  and  vice  versa. 

We  have  therefore  shown  that  the  introduction  of  the  new  class 
of  bo\mdaj-y  conditions  instead  of  the  actual  matching  conditions  at  a 
dielectric  interface  guarantee  a  amall  percentage  error  over  the  entire 
far-field.     The  chief  virtue  of  this  replacement  is  that  it  permits  an 
e:q)licit  approximate  solution  to  otherwise  unsolvahle  problems. 
Consequently  to  gain  assurance  in  the  use  of  the  boundary  condition  it 
would  be  worthirtiile  to  compare  the  result  of  its  use  as  with  the  exact 
solution  of  a  non- trivial  problem.     Such  a  comparison  is  discussed  in 
Section  $, 

5.    Comparison  of  the  Result  of  the  New  Method>  with  an  Sxact  Solution 
for  a  Problem  Involving  Diffraction 

A,     Exactly  Formulated  Problem  and  Solution 

We  consider  the  following  two-dimensional  Wiener-Hopf  problem  for 
which  the  faivfield  in  the  positive  x-direction  can  be  obtained  exactly 
without  need  to  perform  an  explicit  factorization"-  -■ .     A  plane  wave  is 
incident  at  grazing  incidence  along  a  dielectric  interface  (cf.  Figure  7 ), 
For  X  <  0  there  is  a  semi-infinite  metallic  screen  at  the  interface.     If 
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y 


e 


ik.x 


Metallic  Screen  ^^-iz;^ 


Figure  7s     A  plane  wave  launched  on  a  semi-infinite  metallic 
plane  at  the  interface  of  two  dielectrics. 


-  29  - 


we  associate  the  field  u(x,y)  with  the  magnetic  component  H_  then  the 
electromagnetic  boundary  conditions  to  be  satisfied  are 


(1)    S  -  0 
dy 


X  <0,  y  -  0 


and 


(2)    u(x,0+)  -  u(x,0-) 


/-x   1   du(x,0+)   1  au(x.O-) 


X  >  0 

y  -  0 


vrtiere  u(x,y)  is  the  solution  of  the  appropriate  wave  equation 


(U)       (V^  -H  k^)u  -  0 


(5)     (v^  +  k^)u  -  0 


y  >  0 
y  <  0  . 


Ihe  incident  field  is  taken  to  be 


ik^3 


y  >  0    , 


The  boundary  condition  (l)  expresses  the  vanishing  of  tangential  E 
at  a  perfect  conductor.     Equation  (2)  guarantees  the  continuity  of  H_  and 
equation  (3)  implies  the  ccntinuity  of  tangential  E  at  a  dielectric  as  can 
be  seen  from  the  Maxwell  equation 


(U) 


curl  H     «  -itoe  E 
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where  we  assiune  H     =  H     =  0. 
X        y 


To  solve  the  problem  we  note  that  u(x,y)  for  y  >  0  can  be 

r3i 

expressed  in  the  form"-  -' 


"=1=        1       r      A-(v)       i«*i/l4-v' 


where  A~(v)  is  an  unknown  function.     The  branch-cuts  in  the  v-plane  are 
drawn  as  indicated  in  Figure   8  .     Note  that  we  have  asstuned  a  small 
imaginary  part  for  k-   and  kp  for  purposes  of  convenience.     In  a  like 
manner  we  can  express  the  solution  for    y  <  0  as 


z-±y] 


"5 — :? 


J 00               ivx-i  y  kp  -  V    y 
Bllvle^ ^v 
-00            Vkg  -  v'^ 

I 
where  B~(v)  is  an  unknown  function  of  vj  note  the  sign  of  the  radicals 

in  the  exponentials  of  (k)  and  (5).     They  have  been  so  chosen  as  to 

assure  the  outgoing  nature  of  the  solution. 

By  their  construction  the  representations  (5)  and  (6)  satisfy 
the  respective  wave  equations  (3)   and  (li).     In  addition.   If  A"(v)  and 
B"(v)  are  analytic  functions  of  v  for  Im  v  <  0  and  have  the  following 
behavior  in  the  lower  half-plane 


Im 


V  <  0     \ 


|a"(v)|  v-^  |v|  as     |v|^oo  ,  t^>  0 


-6, 


'2 

|b~(v)    \^  |v|  as     |v|— >oo,  Cg  >  0 


However  the  ratio  k^Aj  is  realj   in  the  final  result  these  imaginary 
parts  of  k-    and  k    will  be  neglected. 
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Im  V 


Im  S,  =  + 


Im  $2=  -I- 


ImS  =- 


Im  82=  - 


>  Re 


S,(v)  =  /kf-T^     S,(0)s  +  k, 

=  /k^v^  s(o)= +k^ 


S^iv) 


Figure   fi:  The  complex  v-plane;  the  imaginary  parts  of  S, (v)  and  S„(v)  are 

greater  than  or  equal  to  zero  on  the  real  v-axis.  Two  cuts  in  the 
third  quadrant  are  omitted  in  this  figure,  for  convenience. 
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then  we  can  close  the  respective  contours  for  x  <  0,and  the  boundary 
condition  (1)  will  be  satisfied  by  Jordan's  lemma. 

To  satisfy  the  boundary  condition  (3)  we  note  that  it  is 
sufficient  to  choose 


(7) 


A-(V) 

4 


B"(v) 
,2 


If  Q*(v)  defined 


(8)    -L- .  _£(vL  .  X(v)     .    G*(^) 


by  (8)  is  an  analytic  function  of  v  for  Im  v  >  0  of  algebraic  order  -e-  <  0 
at  infinity  then  the  boundary  condition  (2)  will  also  be  satisfied 
by  Jordan's  lemma. 

It  can  be  shown   [Appendix  k]   that  it  is  possible  to  express 


(9)      P*(v)  P*(-v) 


,-2 


>^    Xi^ 


^(v) 


.  — 1/2 

where  P  (v)  is  an  analytic  function   regular  and  zeroless,  and  0(v     '    ), 
for  Im(v*k;.)  >  0.     We  can  then  rewrite  (8)  as 


(10) 


G^V)    ,      1 
P*(V)        ^-* 


P*(m)      P*(v) 


J^  TO  *  A-(v)  P^-v)   . 
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Owing  to  the  assumed  behavior  of  G  (v)  and  A~(v)  at  infinity  it  follows 
that  each  side  of  (lO)  vanishes  at  infinity.  In  addition  the  left  hand 
side  is  analytic  for  Ira(v+k-)  >  0  and  the  right  hand  side  for 
Im(v-k-)  <  0.  It  follows  that  each  side  is  the  analytic  continuation 
of  the  other.  If  we  invoke  the  Liouville  theorem  we  find  that  each 
side  must  vanish  separately  and  hence 
(U)      A-(v).-^^,ffl(^)   . 

If  we  examine  the  growth  of  0^  (v)  and  P*(v)  we  find 
I  1^    A-(v)  -  0(v^/2j 


so  that  for  y  »  0,  both  u  and  the  normal  derivatives 

s 

'-co 

are  integrable  in  a  vicinity  of  the  origin.  For  y  =  0  the  diffracted  field 
can  be  displayed  explicitly  as 


la.    ./^    P*(ki)o'(v)(v-V 


To  obtain  the  field  for  yO  and  x  »  1  we  note  that  since  o^(v)    is 


♦ 


zeroless  (cf.  Appendix  B)  we  can  deform  the  contour  in  (12)  to  C  as 


-* 


indicated  in  Figure  9  .  The  path  C  is  compared  of  two  portions  t  C^, 
an  Integration  about  the  branch  cut  at  k^,  and  Cg,  another  integration 


about  the  branch  cut  at  kp. 


cf.  Equation  (9)  of  the  present  section,  and  the  last  equation  of  Appendix  A. 
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>Re 


Figure  9:     The  contour  C'*''  is  composed  of  both  C,   and  C 
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For  X  »  1  we  can  obtain  an  aajTnptotic  evaluation  of  (12)  by  the 

[7] 
procedure  of  Van  der  Waerden*-  -■   if  we  discuss  each  integral  separately. 


Consider  first  the  integral  about  C, :     Let  s  -  yv-kT  ,  then 


Kv) 


P*(v) 


Vkl  -V^  (f  (v) 


P*(v) 


^^^ 


i^V^^k^yCp 


admits  a  power  series  expansion  in  s  about  k. 


I(s)  -  Pdt^)    L    +  a^s    +  0 


where 


a^  -  I(k^)  -  1 


and 


,.k. 


-i>M< 


-1/2 


-v^ 


Consequently  the  entire  integrand  of  (12)  ia 


1 


1    *  ^ 


Vv"IS 


9^^   +   0 


(v-k)' 


and  f ollcwing  Van  der  Waerden  we  can  integrate  term  by  term  to  obtain 


(Hi)   u_(x,0)  -  -e 

s 


ik^a 


1  + 


.^rci) 


ni  v^ 


IJ/U 


*     0 
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as  the  contribution  of  the  integration  about  C^  with  the  terms  P(k- ) 
cancelling  each  other. 


The  leading  term  in  an  expansion  about  k^  is  proportional  to 


V^ttt; 


which  contributes  an  amount 


0(:c-V^) 


to  the  asymptotic  development.     Consequently  the  entire  field  for  y  =  0, 
X  »  1  U 


(15)      u(x,0)  -  -    r^        V-P^        ^-7= *  OU'^/^) 

\  *^1         r   k^  -k,  Vnx 


Equation  (15)  gives  the  complete  field  for  jO  and  k^x  »  1  for 
the  problem  of  a  plane  wave  incident  along  a  semi-infinite  metallic  screen 
located  at  a  dielectric  interface.     In  part  B  we  shall  obtain  a  corresponding 
result  when  the  dielectric  interface  is  replaced  by  our  new  boundary  condition. 
It  will  be  seen  that  the  results  compare  very  well. 

B.  Approximately  Refonnulated  Prcblem  and  Solution. 

In  part  A,  formula  (l5)  represents  the  forward  scattering  at  a 
dielectidc  interface  of  a  plane  wave  launched  at  grazing  incidence  on  a 
semi-infinite  metallic  screen.     We  can  obtain  a  corresponding  expression 
to  (1$)  if  we  replace  the  dielectric  interface  by  our  new  boundary  condition. 
The  solutions  for  the  two  problems  can  be  compared  in  the  forward  direction 
(x  »  1,  y»0)  to  test  the  utility  of  the  boundary  condition. 
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The  problem  posed  in  this  part  will  be  identical  with  that 
of  part  A  with  the  exception  that  the  boundary  condition 


("'     (i^#***^^j'"<="^'-<' 


1 

will  replace  the  dielectric  interface  at  x  >  0,  y»0, 

A  and  B  are  the  parameters  introduced  in  Section  1,     In  that  section 
we  indicated  that  there  existed  some  fleocibility  in  the  choice  of  their 
values.     Consequently  we  shall  leave  them  unspecified  except  for  the 
restriction  that  they  be  real  and  satisfy 

A  >  B 
[cf.   (1.8  )]   which  will  imply  the  lack  of  zeros  of  the  expression 


r"?^**-Tr 


k^' 


Once  again  we  introduce  the  construction 
(17)      u(x,y)  .  e^    .  ^   J      -^      e  ^ 


y 
dv 


which  will  satisfy  the  boundary  condition  (1)  if  A."(v)  is  analytic  for 
Im  V  <  0  and, 

(17a)     A"(v)  "^  Ivp*  as  Ivj  -»-oo  Im  v  <  0 


See  Appendix  C. 
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for  some  e  >  0. 

For  comrenience  define 

...  .        ? 

(18)    o^(v)  - 

f,)^?^--tf 

^    2    ' 
k£  -  v^ 

which  can  be  factorized  by  the  well-known  procedure  of  Appendix  A 


(19)        &^{v)  -  P*(v)  p2(-v) 


where  P*(v)  is  an  analytic  function  for  Im(v+k)  >  0  such  that  |P  (v)|  is 
bounded  above  and  below  by  positive  constants  in  its  half  plane  of 
iregnlarity.  It  follows  that  if  G  (v)  defined  by 


(20)   0*(v)  -  |-5|-  +  V^-7  A-(v)  o^(v) 


is  analytic  and  behares  like 


0*(v)  ^  v"S    |v|  ->0  , 


liii  V  >  0 


in  the  upper  half  plane  the  boundary  condition  (l6)  will  be  satisfied (by 
Jordan'  s  lenima]L 

As  usual,  we  shall  now  solve  for  A~(v)  and  G  (v)  by  appeal  to 
Liouville's  theoi^ni.     Rewrite  (20)  as 


G'*'(v)  ^  A-B 

Pt(v)  /^^       ^'h. 


(21) 


A-B 


(v-k^)>^P^(k^) 


>^p2(k^)      VE^  P^(v) 
♦    a"(v)  >^pv  p2(-v) 


At  this  stage  let  k^  have  a  small  positive  imaginary  part  which  will  be 
made  to  vanish  later. 
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Owing  to  the  factor   Vk_-v  of  A~(v)  in  (21)  we  cannot  immediately  conclude 

that  each  side  of  (21)  vanishes.     However  we  have  another  condition  at 

our  disposal:     The  behavior  of  the  solution  in  the  vicinity  of  the  origin. 

2 
We  shall  seek  the  solution  such  that  the  highest  order  derivative  —k  which 

ax 

appears  in  the  bourriary  condition  (16)  be  integrable.     For  this  demand  it 
is  sufficient  that 


|a"(v)|  «^  V  as     |v|  — »'00  ,     Im  V  >  0  , 

where  5  >  1.     A  solution  of  this  type  will  be  obtained  thus  vei-ifying  the 

assumption.     As  a  consequence  it  follows  that  each  side  of  (21)  is  analjrtic 

in  a  half  plane,   shares  a  common  strip  of  analyticity  and  vanishes 

at  infinity  in  its  respective  half  plane.     Invoke  the  Liouville  theorem  and 

obtain 

A-(v)  .         -{^-y L^    p2(.v) 

(v-k^)^  P^(k^)/       ^  * 

or  alternately  by  ((18)  and  (19^  » 

(A-B)  p!   (v) 
A"(v) -j^  *  


(v-k^)^/^o^(v)V3^PA^^^ 
The  integral  representation  of  the  scattered  field  bee  ones 

1^  P;(v)e^^  *i  >^^  7 


(22)      u^(x,y)  -  - 


.00 
A-B  I       "n"  'A'"'"  dv 


2ni  >^  P,(k^)  /^  (^j^)  ^  y^  ,  ,  .^  ^ 


v2 
'1 


Alternatively  we  could  have  required  e  >  1  in  connection  with  the 
equation  (17a). 


-  38  - 


We  can  evaluate  the  far  field  for  x"^;:^  1,  and  y  =  0 
in  the  same  fashion  as  Part  A.  We  obtain 


aa(x,0) 


/+  ivx 

pI(v)    e 
— TTS — 5— 


dv 


A  - 


B  V 


r2     J 


(23) 


-   e 


ik^x 


1.]^  ! 


-TtiA 


A  -  B       ^J^ai 


-3/2 
+   0(x  ) 


Ik^x 
and  the   total  field  e         +  u-(x,0)    is 


(2li.)  u(x 


,0)     =     -  ^     1     e_^ +       Ox 


A-  B 


k  v^ 


for  7=0,   x>$>l. 

Note  that  for  both  the  true  and  the  approximate  problems 
the  leading  term  in  the  scattered  field  cancelled  the  incident 
field*  This  is  of  course  a  conseqence  of  the  fact  that  the  ref- 
lection coefficient  for  grazing  incidence  is  -1  for  both  problems. 
It  is  more  iicportant  to  notice  that  for  both  problems  if  x^^:>>  1, 
the  functional  form  of  the  leading  terras  (1$)  and  (21}.)  agree,  and 
differ  only  by  a  constant  factor 


^1  "H'^l 


1 
k  (A-B) 


*cf.  equation  (15)  of  this  section. 
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or,  since  A_.  ■  l/n,  regardless  of  which  of  the  various  sets  of  (A,B)  we  use^ 


/7^ 


k-,(i-B) 
In 


Consequently  the  percentage  error  between  these  two  coefficients  is 


(25) 


P.  E.  - 


1  - 


vTl) 


¥ 


l-n 


% 


100 


Since  B™,  is  of  the  order  0(l/n-')   [cf .  (1.9)]   for  fitting  at  any  real  angle 
we  see  that  the  error  decreases  with  increasing  n. 

In  line  with  the  remarics  in  Section  1  we  may  observe  that  if 
in  the  far-field  x  »  1,  y*0,  the  most  important  direction  had  been 
defined  as  grazing  incidence  then  one  would  have  chosen  B  so  as  to  match 
the  derivative  of  the  reflection  coefficient  at  grazing  incidence  i.e. 


(26) 


6, 


TM      n 


1  -v5Ti52 


then  there  woitld  have  been  no  error  at  all  in  (2li)  as  compared  with  (15), 


As  we  remarked  earlier  if  we  compare  (2U)  with  (15)  we  notice 
that  the  functional  form  of  the  leading  terms  are  the  same  but  with  differing 
constants  for  various  choices  of  B.  We  can  compute  the  pei^ientage  error 
between  these  leading  coefficients.  In  Graph  9  we  plot  three  curves: 
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(A)  represents  the   "  normal  impedance"    approximation  tt-  ^  +  Au  »  0 
obtained  by  setting  B  =  0  (a  related  condition  customarily  used  to 
represent  almost  perfect  conductors);   (B)  shows  the  error  obtained  by  a 
binomial  expansion  of  the  radical  given  by  (l,ll)j   (C)  is  the  error  if 
B  is  chosen  as  indicated  by  Graph  3  which  is  a  matching  at  Brewster's  angle. 
Notice  that  the  big  improvement  is  from  (A)  to  (b),   an  improvement  ;^ich 
arises  from  the  additional  term  in  the  boundary  condition. 

6,  Conclusion 

To  conclude  Part  I  it  may  be  urged  that  mixed  boundary  conditions 
with  more  than  two  terms  have  a  great  deal  to  offer  in  approximation  techniques 
for  boundary-value  problems.     We  have  introduced  a  technique  by  which  the 
transition  conditions  at  an  interface  between  two  dielectrics  can  be 
conveniently  and  accurately  replaced  by  such  a  single  mixed  boundary  condition. 
The  method  is  useful  when  one  is  interested  in  the  field  in  only  one  of  the 
dielectrics.     We  have  tested  the  accuracy  of  the  results  so  obtainable,   in 
several  ways.     First:  we  have  demonstrated  that  this  is  a  reasonable  and 
well-defined  approach  by  proving  the  uniqueness  and  reciprocity  theorems  for 
a  half  space.     Secondly:  when  a  single  interface  is  present,   so  that  there 
is  no  diffraction,  it  has  been  proved  that  the  maximum  percentage  error, 
over  the  entire  far  field  of  an  arbitrarily  placed  line  source  is  small. 

Additional  confidence  in  the  technique  has  been  demonstrated  by  using 
the  boundary  condition  to   replace  the  transition  conditions  in  a  non-trivial 
problem  which  does  involve  diffraction,     A  comparison  of  the  results  of  the 
approximate  version  and  the  exactly  formulated  problem  has  shown  that  there 
is  excellent  agreement  between  the  two  in  a  certain  respect. 
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Appendix  A 


The  function  ^ 

d(v)  -     ± +  2 


>^^  ^^^4^ 


can  be  factor! zed  aa 

d(v)  =  P*(v)  P*(-v) 

where  P  (v)  is  analytic  and  zeroless  for  Im  v  >  Im  -k,  by  an  application 
of  the  Cauchy  integral  theorem.  Define 

S(V)  8  -^ ^ y^\        


-  Ul  - 


then 

,2       ,2 

o^v)  "   ^        ^        '       S(v)     . 

In  appendix  B  it  is  shown  that  S(v)  and  hence  o^v)  is  zeroless  in  the 
cut  v-plane,  consequently 

/n   S(v) 

is  an  analytic  function  for 

|lin  v|  <  Im  k, 

and  square  integrable  where   |k-|   <  |kp|.     By  the  Cauchy  integral  theorem 

(1)     Jn  S(v)  =  /n  p(v)  -  /nn(v) 

wheire 

+  00  -  ip 

/n  P(v)  =  ^     J  vT^     ^^    S(S^)  dJ^ 

-00  -ip 

and 

+  00  +  ip 


/n  n(v).4     j  J^/nSiK)<^ 


-00  +ip 

where  p  is  any  real  number  0  <  p  <  |k-|.     By  virtue  of  the  evenness  of 
S(v)  we  have  the  relation 


p(v)    »  —7= r 


-  U2  - 

and  by  definition  p(v)  is  analytic  for  Im  v  >  p.  Thus 

S(v)  -  p(v)  p(-v) 
and  if  we  define 


-y^^Tkf  p(v) 


we  have  the  desired  factorization 

d(v)  =  P*(v)  P*(-v)  . 
Since  p(v)  is  0(1)  as  |v|  — >oo  it  follows  that 

P*(v)  -  0(v"^/^) 
as  Ivl— >oo  . 


Appendix  B 
The  expression 

(1)  R(v)-k^yk|T7  +  kl  -Al '  v^ 

can  have  zeros  only  in  a  region  where  the  imaginary  parts  of  the  radicals 
have  different  signs.     In  the  upper  half-plane  Im  v  >  0  this  can  only 
arise  in  the  half  strip  Im  v  >  0, 

Re  k-  <  Re  V  <  Re  kp. 


-  U3  - 


If  we  multiply  (l)  by  R  (v),  where 


(2) 


R*(v)     o    k^Vk^-  v^     -    ^ifkf^^ 


we  find  that  the  roots  of  R(v)-R  (v)  arise  for 


,  2,  2 

k^TT^ 


k''  ITT 

"2 


so  that  V  lies  on  that  oortion  of  the  ray  9  =  &  between  k.  and  the  origin. 

O  '  -L 

But  on  this  segment  a  root  of  R(v)  is  impossible,  since  it  is  to  the  left  of 
both  cuts.     Since  R(v)  is  even  there  are  also  no  zeros  in  the  lower  v-plane. 

Appendix  C 


We  can   show  that  the  expression 


2  (V)   =    ^    ^^^^  .   (A^  -  B^  -^) 
lacks   zeros  in  the  v-plane  if  we  recall   (1.7)  and  (l.ll) 

n 


A    -  :i  -  1  <  1 

TM         kg 


^TM  "  "^TM 


n^  +  1 


-hh  - 


We  have  allowed  k,    and  k^  to  be  complex,   but  in  such  a  manner  that  the 
ratio  k^/kp  is  realj  this  in  turn   implies  that  A       and  B-^^.  are   real. 
Furthermore,   it   is   dear  that  they  are   positive. 

Define  the  expression  Z   (v) 


and  we  shall  show  that    the  zeros  of  the  product  Z(v)»Z   (v)  are  necessarily 
those  of  2*'(v).     For  this  purpose  it   is  convenient  to  introduce  the  new 
variable  4  =  V/4c^  .     We  can  calculate  that  the  product 

Z(v)'Z*(v)   =  Z(k^C)'Z*(k^C)  -  1  -  4^  -  (A^  -  B^  K^) 

is  a  quartic  in  C  and  hence  has  four  zeros.     It  is  easily  shown  that 
these  four  roots  lie  on  the  real  and  imaginary  ^-axes.  To  verify  this 
point  note  that  since  A_,  <  1  the  expression  ZZ     is  positive  at   the 
origin.     As  we  move  along  either  coordinate  axis  in  any  direction  we  find 
that  at  infinity  ZZ     has  changed  sign  and  is  negative.     This  implies  the 
existence  of  two  roots  on  each  of  the  ^-coordinate  axes.     We  have 

^TM  ^  ^TM'    "^"*^®'    '^y  °'^^  choice  of  the  positive  branch  of  the   radical  the 
expression 

has  no  roots  on  the  segment    |P.e  ci   <  1.     Of  course   there  are  also  no  roots 
on  the    segment    JRe  gi    >  1,  since  the  radical  is  imaginary  there.  Hence  Z 


-  h:^  - 


2 
cannot  have  a   zero  on  the    real  ^-axJs.     Since  -?     is  positive  on  the 

imaginary  axis  a   root  of  Z  cannot  exist   there  either.     Hence  all  the 

zeros  of  ZZ     must  be  those   of  Z     and  we  have  proved  the  assertion: 


Z(v)   is  zeroless  in  the  cut  v-plane. 


Appendix  D 

(For  notation,   refer  to 
Section  U) 


Apply  the  wave  equation  to 


Jw^S-v'V" 


and  we  obtain  by  (3.1)  and  (3.2) 


-  U6  - 


i 


V    L 


Up  6(S)   -  Ug   6(R) 


<3K 


Recall  the   fundamental  property  of  the  5-function. 


g(x)  6(x  -  ?)  dx  -  g  U) 


for  any  ranf^e  of  integration  R  including  C«  Hence  as  p 
integral  becomes  Utj(S)  -  Uo(R)  . 

The   integral  about  S,  the  boundary  of  V,  is 


•00  the  volume 


Su. 


Ssr  -^^'  ■^"- 


where  S  is  a  large  semicircle  with  the  x-axis  as  diameter 


and  C,  is  the  circular  portion  of  Sj  we  calculate  that 

/OO  / 


lim 
p  ->oo 


do- 


-00 


dx  ♦lim 

p  ->oo 


dcf 


Consider  the  integral  about  C,    separately,  if  we  assume  that  we  can 


-  U7  - 


interchange   the  limit  and  integration 


lira 
p— >oo 


dvu  bu- 


by  (3.I4) 


=  ik 


(Up  U3   -  U3  u^)  do' 


0     . 


For  the  integral  about  the  x-axis  apply  the  boundary  condition 


(3.3). 


^00 


-00 


/OO 


.^-00 


9^Ug 

Up  (au    *  p  -2-) 
3x 


8  Up 
dx 


dx 


.00  ^2  -2 

/  a    Ug  8    Up 


to  obtain  the  desired  result. 


-  U8  - 
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Attni  John  V.  N.  Granger,  President 

Grumman  Aircraft  Engineering  Corporation 
Bethpage,  Long  Island,  N.  Y. 
Attn:  Mrs.  A.  H.  Gray,  Librarian 

Engineering  Library,  Plant  No.  5 

The  Halllcrafters  Company 

llliOl  West  5th  Avenue 

Chicago  2b,  Illinois 

Attn:  La Verne  LaGloia,  Librarian 

Hoffman  Laboratories,  Inc. 
3761  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Compjany 
Antenna  Department 
Microwave  Laboratory 
Building  12,  Room  2617 
Culver  City,  California 
Attn:  M.  b.  Adcock 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 

Attn:  Dr.  L.C.  Van  Atta,  Associate  Director 
Research  Labs. 

Hycon  Eastern,  Inc. 
75  Cambridge  Parkway 
Cambridge,  Mass. 
Attn:  Mrs.  Lois  Seulowitz 
Technical  Librarian 

International  Business  Machines  Corp. 

Military  Products  Division 
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New  York  33,  New  York 
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Military  Products  Division 
Owego,  Mew  York 
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Attn:  Research  Library 
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Washington  7,  D.  C. 
Attn:  Mr.  Delmer  C.  Ports 
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298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 
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711  Ibth  Street,  N.  W. 
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2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-7?,  Plant  A-1,  Bldg.  63-1 

Purbank,  California 

Attn:  N,  C.  Harnois 

The  Martin  Company 

P.  0.  Box  179 

Denver  1,  Colorado 

Attn:  Mr,  Jack  McCormick 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attn:  Engineering  Library 
Antenna  Design  Group 

Marvland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,  Marvland 
Attn:  Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 

The  W,  L,  Maxscn  Corporation 

U60  West  Tfith  Street 

New  York,  N.  Y. 

Attn:      Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,   St,  Louis  3,   Missouri 

Attn:     R.  D.   Detrich,    Engineering  Library 

McMillan  Laboratory,    Inc. 
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Attn:     Security  Officer.    Document  Rocsn 
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Polytechnic  Research  and  Development 

Company,    Inc. 

202  Tlllary  Street 
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Attn:  Technical  Library 
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Radiation,  Inc. 
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Melbourne,  Florida 
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Radio  Corp.  of  America 
RCA  Laboratories 
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RCA  Laboratories 
David  Sarnoff  Research  Center 
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Attn:  Mr,  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldridge  Corporation 
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Attn:      Engineering  Technical  File 

Page  Communications  Engineers,    Inc. 
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Rantec  Corporation 

Calabasas,   California 

Attn:      Grace  Keener,   Office  Manager 

Raytheon  Manufacturing  Company 
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Attn:      Mr.    Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,    State  Road 
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Attn:      Mr.   Robert  Borts 
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Wayland  Laboratory 
Wayland,    Mass. 
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Republic  Aviation  Corporation 
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Georgia   Institute  of  Technology 
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Harvard  University 
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University  of   Illinois 
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University  of  Illinois 
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Electrical    Engineering 

The  Johns  Hopkins  University 
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Albu-iuerque,    New  Mexico 


Applied  Physics  Laboratory 
The  Johfls Hopkins  University 
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Silver  Spring,    Maryland 
Attn:     Mr.    George  L,   Seielstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  20B-221 
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Attn:   John  H,  Hewitt 

Massachusetts  Institute  of  Technology 
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P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Document  Room  A-229 

University  of  Michigan 
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Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:   J.  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:   Prof,  K.  M,  Siegel 
912  M,  Main  St., 
Ann  Arbor,  Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilantl,  Michigan 
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"niversity  of  Minnesota 

Minneapolis  lb,  Minnesota 

Attn:  Mr.  Robert  H,  Stumm,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E.  Beam 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attn:  Dr,  T.E,  Tice 

Dept,  of  Elec.  Engineering 

The  University  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attni   Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

SS   Johnson  Street 

Brooklyn,  New  York 

Attn:  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
??  Johnson  Street 
Brooklyn.  New  York 
Attn:  Mr.  A.  E.  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec,  Engineering  Research  Laboratory 
P,  '^.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr,  John  R,  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W,  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:   Prof,  G.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R,  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G,  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr,  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:  Prof,  3,  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momingside  Heights 

New  York,  N.  Y, 

Attn:  Dr.  Schlesinger 

McOill  University 

Montreal,  Canada 

Attn:   Prof.  G.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schultz 

The  Pennsylvania  State  ''niversity 
Department  of  Electrical  Engineering 
University  Par'K,  Pennsylvania 

University  of  Pennsylvania 

Institute  of  Cooperative  Research 

3 boo  Walnut  Street 
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Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W,  Cumberland  Avenue 
Knoxville  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 

Seattle,  Washington 

Attn:   Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:   Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 
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200  South  33rd  Street 
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ATTNi  Professor  A,  H.  Waynick,  Director 


Institute  of  Mathematical  Sciences 
25  Waverly  Place 
New  York  3,  New  York 
ATTNi  Librarian 

Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
ATTNi  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.  C. 

ATTN:  Dr.  W.  K.  Saunders 

Applied  Mathematics  and  Statistics  Lab. 

Stanford  University 

Stanford,  Qilifomia 

ATTN:  Dr.  Albert  H.  Bowker 

Departme'^t  of  Physics  and  Astronomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTN:  Dr.  A.  Leitner 

Unlversitv  of  Tennessee 
Knoxvllle,  Tennessee 
ATTN:  Dr.  Fred  A.  Ficken 
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ATTN:  Professor  B.H.  Blssinger 

General  Atomic 
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Department  of  Physics 
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California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
ATTN:  Dr.  A.  Krdelyl 
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Stanford  University 
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ATTN:  Dr.  Harold  Levine 

University  of  Minnesota 
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Department  of  Mathematics 

Stanford  University 
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ATTN:     Dr.  B.  S.   Gourary 
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Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Mass. 
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ATTN:  Dr.  R.  Gallet 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  Qalifomia 

University  of  Minnesota 
The  University  of  Library 
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Department  of  Mathematics 
Unlversitv  of  California 
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ATTN:  Professor  Bernard  Friedman 
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Massachusetts  Institute  of  Technology 

P.  0.  Box  73 
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ATTN:  Dr.  Shou  Chin  Wang,  Room  C-351 

Melpar,  Inc., 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

ATTN:  Mr.  K.  S.  Kelleher,  Section  Head 

Hq,  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTN:  Mr.  Francis  J.  Zucker,  CRRD 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTN:  Dr.  Philip  Newman,  CRRK 


Mr.  N.  C.  Gerson 
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South  Lincoln,  Mass. 

Dr.  Richard  B.  Barrar 
Systems  Development  Corp. 
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Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

UiS  Palisade  Street,  Dobbs  Ferry,  N.  Y. 
ATTN:  Dr.  N.  W.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathanatics 
University  of  Maryland 
College  Park,  Maryland 
ATTN:  Dr.  Elliott  Montroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5,  Mo. 

ATTN:  Professor  J.   Van  Bladel 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street 
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Brandeis  University 
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ATTN:  Ubrary 

General  Electric  Company 
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Stanford  Industrial  Park 
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Electrical  Engineering 
California  Institute  of  Technology 
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ATTN:  DrJSaorges  G.  Weill 

Naval  Research  Laboratory 

Washington  25,  D.  C. 
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Dr.  George  Kear 
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Brooklyn  Polytechnic 

85  Livingston  Street 
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ATTN:  Dr.  Nathan  Marcuvltz 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Jerry  Shmoys 

Department  of  Mathematics 
University  of  New  Mexico 
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ATTN:  Dr.  I.  Kolodner 

Mathematics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Harry  Hochstadt 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 
ATTN:  Dr.  Pullen  Keats 


Dr.  Lester  Kraus 
Ii935  Whitehaven  Way 
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University  of  Minnesota 
Institute  of  Technology 
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Attn:  Dean  Athelston  Spilhaus 

Ohio  State  University 
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Attn:  Prof.  C.  T.  Tai 
Department  of  Electrical  Eng* 
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